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Several families of nonlinear field equations are known to possess space-localized 
singularity-free solutions which describe fields with finite Hermitian norms. This paper 
studies the interaction of such fields ip with given electromagnetic potentials entering 
the field's equations through the covariant derivative forms d^ip — gA^ip . The main 
result is: the motion of a space-localized field, as a whole, is essentially identical to that 
of a point charge moving in the same electromagnetic potentials, if this field is a solution 
either to a nonlinear Schrodinger or to a nonlinear Dirac equation, under some very gen- 
eral and physically meaningful assumptions. Using an alternative method, similar result 
is shown to hold for a very wide class of nonlinear field equations. All of these establish a 
new connection between nonlinear field theory, on one hand, and classical electrodynamics 
and quantum mechanics, on the other. 



1. Introduction 

Since the well known work of Lorentz and Abraham on the spatially extended electron 
several types of space-localized fields (under various names) and their possible role in 
particle/field theory have been studied by many physicists. Among the most prominent 
are: L. de Broglie [1], W. Heisenberg [2], T. D. Lee [3, 4], N. Rosen [5, 6]. . . 

Our interest in the dynamics of space-localized fields, considered as discrete entities, 
is due primarily to the novel link which it establishes between nonlinear field theories, on 
one hand, and classical electrodynamics and quantum mechanics, on the other. In this 
respect the present paper is closely related to an earlier paper [7] by the same author. 

Here we study the global dynamics of space-localized complex- valued multi-component 
fields which are solutions to a large class of nonlinear wave equations. In these equations, 
which are derivable from Lagrangian densities, the coordinates and the time enter explicitly 
only via the covariant derivative forms di^/dx^ -\-iU^{x)ij) , where U^{x) are known real- 
valued functions of the time and space coordinates x M , \i = 0, 1, 2, 3 . Section 3 derives 
the dynamics of a space-localized field, as a single entity, when it is a solution to a nonlinear 
Schrodinger equation modified with the substitutions 

igA^ n = 0,l,2 y 3, (1.1) 



dx^ dx^ 

where g is a real constant and A^ is the 4-potential of a given electromagnetic field. 
Section 4 presents the same derivations in the case when the localized field is a solution 
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to a nonlinear Dirac equation modified with (1.1). The derivations and the results in 
these sections are new. They show that the center of the field's region of localization 
moves as a classical point charge would move in the same electromagnetic field (given by 
Afj,(x) ) provided the latter's variation within the localization region is sufficiently small. 
Section 5 demonstrates that these results are very general by deriving them (using an 
alternative method) from three physically natural assumptions but without assuming that 
the nonlinear field equation is of a particular type. 

N. Rosen [5] (page 98) appears to be the first who has inferred, without giving any 
derivation, a similar result by considering the energy of a space- localized field obeying a 
specific gauge invariant wave equation. Much latter, Bialynicki-Birula and J. Mycielski 
[8] investigating a nonlinear Schrodinger equation with a logarithmic nonlinear term found 
that this equation admits closed-form exact localized solutions, which they called gaussons. 
In the same paper they have demonstrated that if their Logarithmic NLS equation is 
modified with the substitution (1.1) then "...in every electromagnetic field, sufficiently 
small gaussons move like (charged) classical particles." 

Notation and conventions: Here, all complex-valued fields will be denoted with 
the Greek letters ip , (p and r\ , all densities (the integrands of functionals) with script 
capital letters, all 3-vectors with bold letters. Latin indeces will be used for the field's 
components and for the space coordinates. Greek indeces will be used for space-time 
components. The complex conjugate of any field ip will be written as ifj* . The rest of the 
symbols will be defined at their first appearance. The summation convention of repeated 
indexes is assumed as usual. The domain of all space integrals is the entire 3- dimensional 

Q 

Euclidian space R° . The units are selected so that c = 1 . 

2. The assumptions 

Both for clarity and for preciseness the meaning of the term space-localized, as used 
in this paper, will be given by the following: 

Definition 2.1. A scalar, vector or spinor-valued function V = (V'i? • • • > V's) whose 
components Vfe = ^k{x,t) are functions of the coordinates x = (x 1 , x 2 , x 3 ) and possibly 
the time t will be called space-localized, or -field, if its Hermitian norm 



is finite and if b < max x < oo , with b > a constant, for all t ^ . 

All functionals which appear in this paper are assumed to be finite when evaluated 
with a space-localized field. 

It is known that wave equations in 3 space-dimensions which do not depend on x 
explicitly may possess space- localized solutions only if they are nonlinear. Wave packets 
are not localized fields since they dissipate in time, i.e., for them 6 = 0. We will 
assume that in a ^-field there is only a single region of localization, or that if several 
such regions are present they are separated by distances which are much larger than their 




k = 1, . . . , s 



(2.1) 
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characteristic dimensions. Such an assumption will not diminish the generality of the 
following derivations, but it will allow us to focus on the key question: How is the motion 
of the localization region of tp determined by a given external field, i.e., by the explicit 
dependence of the ^ wave equation on x and t ? 

When the -0-field equations are obtained from a Lagrangian density C and certain 
assumptions are met it is possible to derive either from the field equations (Sec. 3 and 
4) or solely from the Lagrangian density C (Sec. 5) the differential equations which the 
coordinates r l = r l (t) of the "center" of localization satisfy. In essence, the first 
assumption is: The external potentials enter the nonlinear field equations exactly as the 
electromagnetic potentials enter the fundamental equations of quantum mechanics. 

Assumption 2.1. The explicit dependence, when present, of the Lagrangian density 

C = C(r, 1>, dr, ty, t, x) = £(V* V, V*, V) (2-2) 
on the coordinates x l and on the time t = x° is only through the forms rj = {rj kfl } '■ 

mix = ^jr + iU^ipk , fe = l,...,s; a* = 0,1,2,3 (2.3) 

where = U^(x,t) are given real- valued functions, if) may be a spinor field, a vector 
field, or just a set of coupled scalar fields, dip stands for the set of all space and time 
derivatives d^ipk = dipk/dx^ . The behavior of and that of the field components ipk 
under rotations in 3 and 4-space are left undetermined in order to maintain generality. 
However, the field equations obtained from (2.2) must be Galilei or Lorentz invariant. 

The above assumption raises the following question: If a nonlinear wave equation 
possesses a space- localized solution, will it possess such a solution after it is modified with 
the substitution (1.1)? This question is answered at the end of this section. 

Assumption 2.2. The functions U^(x,t) are arbitrary except that the magnitudes 
of their variations | AE^| within a localization region are sufficiently small, so that within 
such a region may be treated as being constant. 

A criterion for "when are the variations |AL^| sufficiently small" is given at the end 
of Section 5. 

Assumption 2.3. The Lagrangian density C(ip*,ip,dip*,dip) obtained from (2.2) 
by setting = in (2.3) is gauge type I invariant, that is, the transformation 

Vi = e* £ Vfc, ^ = e- is r k (2-4) 

with real parameter e (which does not depend on x or on t ) preserves the value of C . 

Here, the term "gauge type I transformation" is used as in Goldstein [9] (page 593) to 
distinguish from the gauge transformations of the electromagnetic potentials. 

Assumption 2.3 is a natural extension of a fundamental property of Schrodinger and 
Dirac equations to the nonlinear field equations, considered in this paper. 

When Assumption 2.3 holds it follows from Noether's theorem that there exist a 
conserved density Q and a corresponding flux density T = (JF 1 , JF 2 , J 73 ) associated with 
any C of the form (2.2) and given by 

C = i Lj« * ^ =i U'*-*'X) (2.5) 

V dr] k0 dr] k0 J \ dr) kl dr) ki J 
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where rjko , rjki are defined by (2.3) and tp is a solution of the field equation derived from 
C . In view of the critical role which the expressions (2.5) play here, we will sketch an 
alternative derivation which does not rely on Noether's theorem. 

Let C' = C(ip'*, ip', rj'*, rf) be the Lagrangian obtained from (2.2) by transforming 
the V'-field according to (2.4). Then, the invariance condition is 



dC 



de 



, dC ,„ dC dC „ dC , , 

... = **' W k ~ M + ~ m »s^ = • (2 ' 6) 



Adding to this equation the Euler-Lagrange equation 

dC . TT dC d dC 

multiplied by —iipk and the complex conjugate of (2.7) multiplied by produces 
immediately the conservation law 

dQ 

fL±+y.F=0 (2.8) 

and the expressions (2.5) for which (2.8) holds. 

In what follows we will need a variable r = (r 1 , r 2 , r 3 ) G IR 3 which gives the position 
of the localized ^-field as a whole. This variable is defined, just as the position of any 
distribution, with the functionals 



r l = ^ j Qx l d 3 x , with Q = J Qd 3 x , i = 1,2,3 (2.9) 



~ Q 

R 3 R 3 



where Q is given by (2.5) and Q is a constant of the motion according to (2.8). The 
values of the functionals r l are functions of t only and may be regarded as the coordinates 
of the "center" of the localization region. No probabilistic interpretation will be imposed 
on the density Q , or on any other density, since here such would be entirely unnecessary. 
The integral in (2.9) for r l is normalized by dividing it with Q since the -0-field, being a 
solution of a nonlinear equation in general, will not remain a solution if it is normalized. 
The curve r = r(t) will be called trajectory of the localization region. 

If a "free" wave equation obtained from the Lagrangian (2.2) with = possesses 
a localized solution, when are arbitrary functions no localized solution will exist, in 
general. However, if are strictly constant, but otherwise arbitrary, within the space 
occupied by the •0-field then a localized solution will certainly exist. To see this, apply to 
the Lagrangian (2.2) the transformation ip' k = exp(z#) ipk m which 9 = U^x^ (not a 
Gauge type I transformation), where all = const . Then, observing that the forms 
rjk^ (2.3) when written in terms of ip' k become 

9ipk ,. TT i -id d ^k 
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we apply Assumption 2.3 (for which the dependence of 9 on x and t is now irrelevant) 
to obtain the equality 

C{^*,^,rj*,ri) =£(e l V* e"*V, e id d^'*, e~ ie dtf) = C(ip'*, ip' , dip'*, dip') . 

This shows that if we set = in the equations for ip (but not in # ) we get the 
equations for tp' . Recognizing that the Hermitian norm of ip is not affected by the 
transformation ip' k = exp(i6>) ipu the above claim follows. 

Now, continuity considerations lead us to conclude that even when are not strictly 
constant the localized solution will still exist, as long as the variation's magnitudes | AL^| 
within the region of localization remain smaller than a certain critical value k\ . This is 
one reason for having Assumption 2.2. The second reason is that a Lagrangian function 
for the coordinates r l {t) can be defined, as we will see in Section 5, only when |AL^| is 
smaller than some other limit value k^ ■ It may be argued that k\ = k<i , but for our 
purposes this is not needed. 



3. The localized field obeys a nonlinear Schrodinger equation 

Consider the family of nonlinear Schrodinger (NLS) equations 

i?l ¥ = "^ v ^^ W) (3 - 1} 

for the scalar- valued field ip . The various constants which appear in Schrodinger equation 
proper have been retained in (3.1), so that a comparison with certain results from QM 
can be made most conveniently. It is known that there are functions G = G(ip*ip) with 
which (3.1) possesses space-localized solutions. For the existence conditions and the proofs 
see Berestycki and Lions [12]. 

The Lagrangian density for the NLS equations (3.1) is 

where G(p) = dQ(p) / dp . In the presence of an electromagnetic field the substitution 
(1.1) with g = q/h is inserted in (3.2), according to Assumption 2.1. Then, equation 
(3.1) becomes 



( — V-V + g$ + G(V>*V)W , V = (ihV + qA) (3.3) 

which will also have localized solutions if Assumption 3.2 is met (see the discussion at the 
end of Section 2). The first of expressions (2.5) gives Q = ip*ip . Thus, the coordinates 
r % of the center of the localization region, according to (2.9), are 

r l = -1 J ^> x i d 3 x , N = J d s x . (3.4) 

R 3 R 3 
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With the use of (3.3) the velocity v = (v 1 , v 2 7 v 3 ) of the localization region is 
dr Iff,* dip dtp* , \ )3 

v= Tt = N ^ x m + -dt x V dx 



[V-V,x] + [G,x]^d 3 x = Jij*Vijd 3 x (3.5) 



TiN 

since the commutator [V • V, x\ gives 

[V • V, x] = [V j V j , x] = V j [V j , x] + [V j , x] V j = — V 

and since the function G(ip*ip) certainly commutes with x . 

According to (3.5) V is the velocity operator for a NLS equation and it is the same 
as the velocity operator in QM. We use again (3.3) to find the second time derivative of r 

d 2 r dv 1 [{ dip dtp* dV \ 3 

d^ = ti = Nj r V ^ + ^T V ^ ~WV dx (3 ' 6) 

For any space-localized function ip the integral containing the nonlinear term is identically 
zero, as seen from the following 

JiP*[V,G]ipd 3 x= 1 -^- J(ij*GVip-ip*V(GiP))d 3 x= 1 -^- jGV(*p*ij)d 3 x 

- S / nzk?™ d3x = £ / v w> A = ■ (3 - 7) 

One short and one long calculations, both routine in QM, yield for the remaining two 
commutators: 

[V, $1 = - — V$ (3.8) 
L J m 

[V,V-V] =V j [V,V j ] + [V,V j ]V j = 4(VxB-BxV) . (3.9) 



When (3.7), (3.8) and (3.9) are inserted into (3.6) and E = -V$ - dA/dt is taken into 
account the result is 



d 2 r q 



y^*(E + ^(VxB-Bx vf)ipd 3 x . (3.10) 



It is most remarkable that the nonlinear term G(ip*ip) of the NLS equation does not 
enter in the above equation. Hence, we would have gotten exactly the same result had 
we taken ip to be a solution of the Schrodinger equation proper ( G = ). However, the 
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derivation of equation (3.10) is valid only if ip is a space- localized field. But, in quantum 
mechanics the only space-localized fields are the bound states in some potential well, and 
so, the "center" coordinates r (3.4) of any such bound state are fixed by the potential 
well location. The wave packets of QM are not localized fields, as explained in Section 2. 
Thus, the question "what is the motion of a localized field as a whole?" is not meaningful 
in quantum mechanics. 

For the space- localized solutions of a NLS equation the above question is certainly 
meaningful. Its answer would follow immediately from (3.10) had we assumed that the 
variations of the fields E , B , instead of those of the potentials $ , A , are negligible 
within the localization region. (Why Assumption 2.2 was chosen to restrict the potentials 
and not the fields intensities will become clear in Section 5.) The equations r = r(t) for 
the motion of the localization center do follow from equation (3.10) and Assumption 2.2 
but the reasoning is not entirely straight forwarded: We change the integration variable x 
to x' according to x = x' + r , where r is given by (3.4). Then, the first term in (3.10) 
becomes (with the abbreviations <&(cc) = <&(cc,£) and A(x) = A(a?,t) ) 

J Tp*ip E(x) d 3 x = -J tfj*tfj(v x ^(x) + Q* ) d 3 x 
= -V r j ip*ip ®(x' + r) d 3 x' - j i>*i> ^A(a/ + r) d 3 x' 

= - V r $(r) J d 3 x' - ^A(x' + r) d 3 x' + J A(x' + r) * f d 3 x' 

= -V r $(r) J ^^d 3 x - J ip*ipd 3 x = E(r) J ip*ipd 3 x = NE(r) (3.11) 

since within the localization region $>(x' + r) = <fr(r) and A(x' + r)=A(r) according 
to Assumption 2.2. Similarly, the term containing B in (3.10) becomes 

^ |f (VxB-BxV)^ 3 x = -B(r)x Jtjj*Vijd 3 x = NvxB(r) (3.12) 

where v = dr/dt , given by (3.5), is the velocity of the localization center. Inserting 
(3.11) and (3.12) into (3.10) produces the desired result 

m-L = / f^ 3 i-Bx U*Vijd 3 x} = q(E(r) + exB(r)) (3.13) 

which is the equation for the localization center coordinates r(t) . This equation is identical 
with the equation for the motion of a classical point charge q with mass m in given 
electric field E(r) and magnetic field B(r) (see Jackson [10]). 



4. The localized field obeys a nonlinear Dirac equation 

Now let ip = (ipi , . . . , ip4) be a 4-spinor field which satisfies a nonlinear Dirac (NLD) 
equation of the form 

^ft7 M ^ - (Mi + 7 5 M 2 )^ = (4.1) 
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where Mi = Mi(si, s 2 ) and M 2 = M 2 (si, S2) are scalar real- valued functions of the 
scalar and pseudo-scalar quantities 

si = ipip and s 2 = ip^ip ( V> = 7/^7° ) . (4.2) 

As explained in Section 2, all wave equations considered in the present paper must be 
derivable from a Lagrangian density. If one takes for the above two functions: 

Mi= dM( Sl ,s 2 ) ^ M2 = dM{s u s 2 ) 
as\ ' 0S2 

then, the real-valued Lagrangian density from which (4.1) is obtained has the form 

where M. = JA(si,S2) is some scalar real- valued function of si and s 2 and nothing 
else. It is known that, with certain M , (4.1) possesses space-localized solutions: Cazenave 
and Vazquez [13], Soler [14]. Additional references may be found in both of this works. 

When external EM potentials $ = — Aq , A = (A\, A 2 , As) are present the equation 
(4.1) is modified with the substitution (1.1), according to Assumption 2.1. Space-localized 
solutions will exist if the original equation (4.1) possesses such solutions and if Assumption 
2.2 is met (by the argument discussed at the end of Section 2). Equation (4.1) thus modified 
and written as a Hamiltonian evolution equation becomes 

ih^- = -a- (ihV + qA)ip + q$ip + Mip = (a -11 + q§ + M)ip (4.5) 

where 11 = — i%V — qA , a = (a 1 , a 2 , a 3 ) with a 1 = 7°7 l , and 

M = 7 (Mi +7 5 M 2 ) . (4.6) 

The various constants in equation (4.5) are written just as they appear in its linear version 
so that a comparison with certain results in Dirac's theory can be made most conveniently. 
We do the substitution (1.1) in the Lagrangian (4.4) and insert the result into the expres- 
sions (2.5). Then it follows that the norm iV , given by (2.1), is a constant of the motion 
as it is in the linear case. 

The position of the localization center r = (r 1 , r 2 , r 3 ) is defined again by 

T% = iv 7 / ^kx*d 3 x = iP^x'd^x . (4.7) 

R 3 R 3 

Then, the velocity of the localization region v = (v 1 , v 2 , v s ) , with the use of (4.5), is 

^4/(^f*>) A 

= -i/*V,„.V]^=i/*taV^ (4.8) 



8 



Complex Space-localized Fields 



Bodurov 



which is the same as the expectation value of the velocity operator a in Dirac's theory 
because the function M , as a multiplication operator, obviously commutes with x % . 

In the non-relativistic case the "mechanical" momentum of a space-localized field is 
obtained, of course, from the velocity functional (3.5) by multiplying it with the mass m . 
In Dirac's theory one takes for the "mechanical" momentum exactly the same functional, 
except in it ifj is a spinor- valued field. Here we do the same, i.e., we take 

p= J ^(ihV + qA)^d 3 x = ^U^d^x (4.9) 

to be the "mechanical" momentum of the localized field, where II = — ihV — qA is 
the corresponding operator in both the relativistic and non-relativistic QM, Messiah [15]. 
(The term "mechanical" momentum is used as in references 9 and 15.) Now, instead of 
calculating dv/dt which was done in Section 3, here we calculate dp/dt as follows 

= ^/V([a-n,n] +«[*,n] + [M,n\+i q h^-)^d 3 x . (4.10) 

For the first and the second commutators we give only the final results since they are well 
known in QM, i.e., 

[a -11,11] = 0? [n^n] = -iqhax (V x A) = -iqhax B (4.11) 
[$,n] = zftV$ . (4.12) 

Although the calculation involving the nonlinear term ip^ [M, Tl\ip is not complicated it 
is presented below in some detail because it is critical for the present derivation. Indeed, 
the following will show that the integral of this term over all space is identically zero: 
Substituting the expression for M from (4.6) and taking in consideration (4.3) 



^[M,n]i(;d 3 x = -ih J (^MVf + {Vip ] )MiP) d s x 
= t/ ((WV^ + ^VV^Mi + ((VV' t )7°7 5 ^ + V' t 7 7 5 VV')M 2 ) d £ 
= -J (MxVtyVV) +M 2 V(^7°7 5 V)) d 3 x = -J (§^V*i + §^Vs 2 )d 3 x 



l3 x 



J\/M(s u s 2 )d 3 x = (4.13) 



one arrives at the desired result, where as before si = and s 2 = V'W 7 5 V' • 

Inserting (4.11), (4.12) and (4.13) into (4.10) and noticing that E = -V$ - dA/dt we 
finds that 

^ = ± [^(E + axB)if;d 3 x (4.14) 

(Jib i V J 
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which is valid for the nonlinear Dirac equations (4.5), as well as for the Dirac equation 
proper. This is so because the integral (4.13), involving the nonlinearity M , is zero. 

Assumption 2.2 was not needed for the derivation of (4.14). Now, if we assume that 
the variations of $ and A are sufficiently small within the localization region, i.e., if 
Assumption 2.2 holds, then the electromagnetic fields E , B may be treated as being 
constant in that region, as it was shown in Section 3. Hence, (4.14) becomes 



dp q 
~dt~N 



(E(r) J ^if; d s x - B(r) x J tp d s x^j = q (E(r) + v x B(r)) (4.15) 



where v = dr/dt, given by (4.8), is the velocity of the localization region "center". 
Equation (4.15) describes the motion of the localization region r = r(t) in complete 
agreement with the result (3.13) of Section 3. 

Again, one recognizes this to be precisely the equation of motion for a point charge q 
but this time in the relativistic case (Goldstein [9] Sec. 7.8, Jackson [10] Sec. 12.1). 



5. The Lagrangian function for the motion of a localized field 

The results in the previous two sections are valid for entire families of NLS and NLD 
equations which possess space- localized solutions. Thus, the question arises: Is it possible 
to derive similar results without assuming a particular type for the field equations? This 
section will show that it is possible, but the method has to be different since one cannot 
calculate explicitely the time-derivatives of the functionals r nor of v without knowing the 
field equations. Instead, we will derive from the assumptions in Section 2 the general form 
of the Lagrangian function for the motion of a space-localized field, i.e., for the coordinates 
r = (ri,r2,rz) of the localization center. Then, several important conclusions will be 
drawn from the form of this Lagrangian function. 

When Assumption 2.2 holds, i.e. the variations |AL^| of the given potentials 
Ufj, are sufficiently small within the localization region, the value of the space integral of 
the Lagrangian density (2.2) may be approximated by replacing in the latter the actual 
potentials U^(x,t) with their values U^r, t) at the localization "center" r. This is so 
because C — > outside the localization region. Then, at a fixed time t and for a fixed 
space-localized tp(x,t) the Lagrangian functional 

L(U(r,t))= j C{i}*,^,d^* -iU(r,t)if>*,dil> + iU(r,t)il>)d 3 x (5.1) 

R 3 

becomes a function of . Alternatively, one can regard (5.1) as a functional of ip and ip* 
in which enter as parameters . At the end of this section it is shown that the error in 
this approximation can be made arbitrarily small by choosing potentials whose variations 
within the localization region are sufficiently small, i.e., |AL/q| + • |AL/j| <C \L/Q\ . 

The function L(U) defined by (5.1) is different, of course, for different ip . Our goal is 
to derive a relationship between L and which remains valid during the entire motion 
of a particular ip as long as ip remains space-localized. To do this we differentiate (5.1) 
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with respect to . Because of Assumption 2.2 it is permissible to differentiate under 
the integral sign which produces 

dL f dC ,„ . f ( , dC , * dC 



This is yet another functional of ip in which enter as parameters. However, if tfj is 
a space-localized solution of the field equations (2.7), derived from the same Lagrangian 
density J0,(ip*, ifi, rj*, rf) , and are the values of the potentials at the region of localiza- 
tion, then, the integrands in (5.2) are precisely the conserved density —Q and the three 
components of the flux density —T % as given by (2.5). 

Thus, when (5.2) is evaluated on a solution for fx = we have that 



dL 



Qd 3 x = -Q (5.3) 



holds during (and only during) the motion, where the quantity Q is a constant during 
the same motion according to (2.8). Of course, the value of Q is determined by the initial 
conditions of the particular solution and so Q will be different for different solutions V • 
Again, when (5.2) is evaluated on a solution but now for fx = i = 1, 2, 3 we have 

Next, the velocity v = (y 1 , v 2 , v 3 ) of the localization region is expressed in terms of Q 
and T % as follows: We differentiate the position functionals r % in (2.9) with respect to t 

i dr l 1 fdQ lj3 1 [dP lj3 1 f . 3 

v l = — = — / — x l d 6 x = - — / r- x l d 6 x = — / T l d s x (5.5) 

dt Q J dt Q J dx 3 QJ 1 ' 

then we use equation (2.8) and integrate by parts. The last result together with (5.4) 
yields a relation which holds for the entire motion 

— = -Qv\ i = 1,2,3. (5.6) 

Let us now consider the motion of the localization region of a solution ip = ip(x,t) , 
assumed to be known, of some nonlinear field equation. We can find the exact trajectory of 
the localization region r = r(t) by inserting this solution into the functionals (2.9). At 
any point r(t) on this trajectory the quantities Uq = Uo(r,t) and Ui = Ui(r,t) which 
appear in (5.3) and (5.6) are well defined and known, since U^{x,t) are given functions. 
Furthermore, the velocity v(t) = dr(t)/dt is well defined and known, since r(t) is known. 
We want to find a function L = L(U,v) of the variables Uq , Ui, v k (i,k = 1,2,3) 
such that the equations (5.3) and (5.6) are satisfied during the motion, i.e. at each point 
of the trajectory r = r(t) . The unique solution of this system of 4 equations, obtained 
by four consecutive (trivial) integrations with respect to Uq , U\ , U2 and U3 , is 

L = -QiUo + Uv 1 ) +L(0,v) (5.7) 
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where L(0,v) depends on v l but not on and may be regarded as the integration 
"constant". One can verify the expression (5.7) by inserting it in (5.3) and (5.6). 

It remains to convince ourselves that L given by (5.7) is the Lagrangian function for 
the trajectory r = r(t) . This is done by restating (5.7) as 

/ £(V>* V, V*, V) d 3 x = L(0, v)-Q (U + U t v l ) (5.8) 
Jn 3 

and integrating both sides with respect to time from t\ to ti . Thus, within the error of 
the approximation the equation 

r I C(iP*,iP,r l *,r l )d 3 xdt = [ \l(0, v) - Q(U + U t v*)) dt (5.9) 
Jtx Jn 3 Jt 1 

holds on space-localized solutions ip(x, t) . Moreover, this equation holds also for all space- 
localized functions ip(x, t) which satisfy the conservation law (2.8) and Assumption 2.2, 
because for such functions (5.8) holds as already shown. Pertinent remarks: 

(a) ip will be called the varied field if it is sufficiently close in norm to a solution ip . 

(b) ip is restricted to be a space-localized field since otherwise the Lagrangian functional 
J Cd 3 x when evaluated with ip will be divergent, in general, and L will be meaningless. 

(c) The varied field ip must be such that it represents a particle whose trajectory is varied, 
but not its properties like mass, charge, etc. Hence, the variations are restricted so that ip 
satisfies the conservation law (2.8) and Assumption 2.2. Consequently, Q in (2.9) remains 
a constant of the motion (which is not varied). 

(d) The functions Uo(x,t) and Ui(x,t) are not varied — they are given. 

According to the calculus of variations, ip(x,t) is varied so that its values ip(x,ti) 
at ti and ip(x,t2) at t% remain fixed, i.e. ip(x, t±) = tp(x, t\) and ip(x, £2) = ip(x, £2) • 
Inserting the varied ip into the expressions (2.5) and the result - - into (2.9) one finds the 
varied trajectory r = r(t) . Thus, by varying the localized field ip we vary the trajectory. 
Moreover, since the functions ip(x, t\) and ip(x, £2) are fixed it follows from (2.9) that the 
initial point f l (ti) and the final point f 1 ^) of the varied trajectory are also fixed, i.e. 
r l (ti) = r l (ti) and P(t2) = ^(£2) , just as required by the calculus of variations. 

Having r l (t) one can calculate the localization center velocity v l (t) and the values of 
Uo(f, t) , Ui(r, t) on the varied trajectory. Finally, one can calculate the value of L(0, v) 
on the same trajectory by setting Uq = and Ui = in (5.1), i.e. from 

L(0, v) = J C(ip*, ip, dip*, dip) d 3 x . (5.10) 

Putting all the above together we see that the expression 

L(0, v)-Q {U (f, t) + PUtf, t)) (5.11) 

has a well defined value at each point of the varied trajectory. Consequently, the value of 
the time integral of (5.11) from t\ to £2 is well defined and equation (5.9) holds with any 
varied field ip under the above stated restrictions. 
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Now, consider a ip for which the left-side integral in (5.9) is at a minimum (maximum). 
According to the calculus of variations this ip is a solution of the corresponding Euler- 
Lagrange equations. Then, as long as equation (5.9) holds, the right-side integral is also 
at a minimum (maximum) which is achieved for the "true" trajectory r(t) . This shows 
that the right-side integral in (5.9) is the action integral for the trajectory r(t) of the 
localization center and that L given by (5.7) is its Lagrangian function. End of proof. 

If we set Up = gA^ in the Lagrangian function (5.7) we see that the first term in it 

-gQiAo + v'Ai) (5.12) 

is the part which accounts for the interaction between the space-localized field , as a whole, 
and the electromagnetic field . This part is identical with the interaction term in the 
Lagrangian function for a point charge q = gQ , both in the classical and in the relativistic 
cases. For the Lagrangian of a point charge see, for example, Goldstein [9] (Section 1.5) 
or Jackson [10] (Section 12.1). The second part in (5.7), clearly, stands for the Lagrangian 
function for a "free" localized field, as a whole, and it is a function of v only. To find it set 
= in the field equation (now assumed to be known) and obtain the corresponding 
stationary localized solution. Apply to it a Galilei or Lorentz transformation to obtain a 
solution whose region of localization is moving with some constant velocity v . Evaluate 
the Lagrangian functional (5.1) with this transformed solution and with = . Since 
v l are the transformation parameters, the result will be a function of v , as claimed, and 
nothing else. We summarize the first principal result which follows from (5.7): 

If Assumptions 2.1, 2.2 and 2.3 are met the interaction term in the Lagrangian 
function (5.7) for the motion of a space-localized ifj -field, as a whole, is identical to the 
interaction term in the Lagrangian function for a point charge proportional to Q moving 
in an electromagnetic field proportional to , regardless of the type of wave equation 
which ifj satisfies. 

For the second result which follows from (5.7) we observe that the above statement is 
valid, in particular, when ifj satisfies a Hamiltonian evolution equation of the form 



1 = T7T wlth 

dt 5ip* k 



H= J Hty*, ip, rf, V) d 3 x (5.13) 



where H is the Hamiltonian density, rjk^ are given by (2.3) and SH/Sift^ is the 
variational derivative of the Hamiltonian functional H with respect to Vfc • F° r details 
see Bodurov [7]. Indeed, it is easy to see that (5.13) is the Euler-Lagrange equation derived 
from the Lagrangian density (which satisfies Assumptions 2.1, 2.2 and 2.3 if H does): 

When C (5.14) is inserted into (2.5) it yields 

Q = ^ k and Q = (V>, V) = N . (5.15) 

The form (5.13) contains, as special cases, Schrodinger and Dirac equations together 
with their nonlinear versions. In these equations the constant g has the value g = q/K . 
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On the other hand, the interaction term (5.12) and expressions (5.15) tell us that the 
electric charge is q = gN . Hence, we have q/Ti = g = q/N from which follows that 
h = N . The dimensionality of N (ip is not normalized) is certainly right, i.e., [action] 
as seen from (5.13) by noticing that the dimensionality of H is [H] = [energy] /[volume] 

[N] = [\ip\ 2 ] [volume] = [H] [time] [volume] = [energy] [time] = [action] . 

Thus, our second main result is: 

If a space-localized ip -field is a solution of a nonlinear complex Hamiltonian evolution 
equation with the form (5.13) whose Hamiltonian functional is gauge type I invariant, then 
in certain relationships the Hermitian norm (ip, tp) appears in the place in which Planck's 
constant fi appears in the corresponding relationships of quantum mechanics. 

Two arguments independent from the one above and pertaining to different relation- 
ships but leading to the same conclusion above, were presented in a paper by the same 
author (Bodurov [11], Sections 3 and 4). A year latter, one of these was discussed with 
further details in Bodurov [7] (Section 5). 

The criterion for the validity of Assumption 2.2, i.e., for the maximum rate of change 
of U at which equation (5.7) still holds is derived as follows. Let a be an estimate of the 
radius of localization given by the functional 

a 2 = ^Yl [ Q(x l ~ r 1 ) 2 d 3 x (5.16) 

where r % are the coordinates of the "center" of localization given by (2.9). Denote by 
t^max and t^min the values of U , within a sphere with a radius a and a center r , for 
which the magnitude \L(U mgx ) — L(U m i n )\ is maximal. L(U max ) and L(U m i n ) are the 
values of the Lagrangian functional (5.1) when it is evaluated with the constant U = ?7 max 
and U = U m i n . Then, the sought criterion for AU^ = (L^) max — (L^) min is 

|L(C/ max ) -L(U miQ )\ < \L(U miQ )\ . (5.17) 

(We assume, without loss of physical generality, that are monotonic functions of the 
coordinates within the localization region.) With the use of (5.3) and (5.6) we can write 



AU^ = -Q(AU + v i AU i ) 



The last equation together with (5.17) and (5.12) produces \AU \ + \v l \ ■ \AUi\ 
\L(0,v)/Q\, where \AU \ and |AL^| are the variation's magnitudes of the respective 
potentials over the localization region, v l are its velocity components and L(0, v) is the 
value of the Lagrangian for the "free" -0-field. The above inequality is equivalent to 



\AU \ < 



L(0,v) 



Q 



and K| • \AUi\ < 



L(0,v) 



Q 



1,2,3 (5.18) 
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that is, the interaction terms in the Lagrangian (5.7) must be much smaller than the "free" 
•0-field Lagrangian L(0, v ) — a physically meaningful conclusion. 

6. Concluding discussion 

It was demonstrated in this paper that if the interaction of a nonlinear space- localized 
V'-field with external electromagnetic potentials is defined according to the substitution 
(1.1) and Assumptions 2.1, 2.2 and 2.3 are met then the motion of the ifi-Reld, as a whole, 
is identical to that of a classical point charge. Moreover, this identity is largely independent 
of the particular -0-field equation, as long as the latter is derivable from a gauge I invariant 
Lagrangian. 

The results (3.10) and (4.14) may be regarded as generalizations of Ehrenfest theorem 
showing that this theorem is valid not only for the linear equations of quantum mechanics 
but also for the families of NLS and NLD nonlinear field equations which possess space- 
localized solutions. 

One may also say that according to the results (3.13), (4.15) and (5.12) an external 
electromagnetic field given by a 4-potential, under the Assumptions 2.1, 2.2 and 2.3, acts 
with the Lorentz "force" gQ(E + wxB) on a space-localized ip-Held, whose equivalent 
electric charge is q = gQ . 

As already discussed, the global description of the time evolution of a ip-Reld by the 
trajectory of its center of localization r = r(t) is valid when the variations | AL^| within 
the localization region are sufficiently small. It is worth noticing that this correlates very 
well with a number of experimental evidences that the concept of elementary particle, as 
an individual entity, is meaningful only when the intensities of the external fields acting 
on this particle are within certain bounds. 

We would like to look on the results summarized above as new evidence that it is not 
only possible but it is also very desirable to represent elementary charges with localized 
•0-fields. For, these results point to a unique opportunity to reconcile the discrete aspect 
with the wave aspect of elementary particles/charges. 
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